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Abstract 

A dynamical correspondence is established between positively curved, isotropic, 
perfect fluid cosmologies and quasi-two-dimensional, harmonically trapped Bose- 
Einstein condensates by mapping the equations of motion for both systems onto 
the one-dimensional Ermakov system. Parameters that characterize the physi- 
cal properties of the condensate wavepacket, such as its width, momentum and 
energy, may be identified with the scale factor, Hubble expansion parameter 
and energy density of the universe, respectively. Different forms of cosmic mat- 
ter correspond to different choices for the time-dependent trapping frequency 
of the condensate. The trapping frequency that mimics a radiation-dominated 
universe is determined. 
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In recent years, analogies between various condensed matter systems and different 
branches of gravitational physics have been developed. (For reviews, see, e.g., Refs. 

For example, the propagation of acoustic waves in an irrotational, inviscid, 
barotropic fluid is formally equivalent to that of a massless scalar field on a curved, 
Lorentzian spacetime This implies that quantum field theory in curved space can 
be modeled in terms of the quantization of the sound wave (phonon field). Further- 
more, it is possible to model a black hole acoustically in terms of supersonic fluid flow 
and, in principle, quantum effects associated with black hole event horizons may then 
be studied within the context of condensed matter configurations [3J EJ |5] . A related 
phenomenon is that of cosmological particle production arising in a time-dependent 
gravitational field and the description of such a process in a dynamical Bose-Einstein 
condensate has recently been investigated [UJ |Hl E] . 

Identifying a link in this way between gravitational and non-gravitational physi- 
cal systems is of great importance. To date, attention has focused primarily on the 
kinematical properties of general relativity, such as the existence of event horizons. 
The purpose of this work is to note that a dynamical correspondence may also be 
established between isotropic, four-dimensional cosmological models and harmon- 
ically trapped, quasi-two-dimensional Bose-Einstein condensates. The correspon- 
dence arises because the equations of motion for both systems can be mapped onto 
the one-dimensional Ermakov system ^UJ ■ Ermakov systems arise in a variety 
of physical phenomena ^21 El El and the one-dimensional case corresponds to a 
second-order, non-linear, ordinary differential equation (ODE) |TTH lllj. 

We begin by considering the condensed matter system. A Bose-Einstein conden- 
sate is the ground state of a collection of interacting bosons trapped by an external 
potential. In the limit where the number of atoms is sufficiently large and the atomic 
interactions are sufficiently weak, the mean-field approximation may be invoked, 
where the effect felt by a particular atom due to the ensemble is approximated by the 
mean action of the entire fluid on the particle. (For a review, see, e.g., Ref. [Hi])- In 
this case, the macroscopic wavefunction for the condensate, is determined by the 
Gross-Pitaevskii equation ^H] 

lh ^Z = + V ( r > *)* + (!) 

ot 2m 

where m is the mass of the atoms in the condensate, V(r, t) represents the trapping 
potential, and g parametrizes the strength of the atomic interactions. 

Eq. (P) is formally equivalent to a three-dimensional, non-linear Schrodinger 
equation and, in general, it is non-integrable. However, in two dimensions it is 
possible to determine the dynamics of the wavefunction by employing the 'moment 
method' ^7J E5] . In this approach, integral relations are constructed directly from the 
wavefunction without solving the Schrodinger equation explicitly and the evolution 
of these physical quantities then parametrizes the dynamics of the wavepacket. 

In view of this, we consider a two-dimensional, cylindrically symmetric Bose- 
Einstein condensate in a parabolic trapping potential V(r,t) = muj 2 r 2 /2 with a 
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time-dependent frequency, u = u(t). Lower-dimensional Bose gases are currently 
of interest both from the theoretical [213 HH E3 122] and observational perspectives 
and two-dimensional condensates have recently been observed experimentally [2H] • A 
Bose-Einstein condensate becomes effectively two-dimensional when the excitations 
in a given direction are frozen and this occurs when the confining frequency in that 
direction is significantly greater than the energy for the thermal fluctuation. 

We now briefly summarize the moment method [17] ITS] ITU] . For the case of a 
constant interaction, g, the four integral parameters [T7] IT9] 



h{t) = 
h(t) 



h(t) 
hit) = 

i I d 2 xr 

d 2 x 



d 2 x\V\ 2 
d 2 xr 2 \^\ 2 



\j/ — 

dr dr 
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are defined, where it is assumed that ^ = ^(r, t) and we have specified H = m = 1 
without loss of generality. Eqs. (JSjH© admit physical interpretations in terms of the 
norm, width, radial momentum and energy of the wavepacket, respectively 

By differentiating Eqs. in turn and substituting in Eq. (JJJ), it can be 

shown that these integral parameters satisfy the set of coupled, first -order, ODEs 

m 



dh 
~dt 
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It may now be verified by direct differentiation that the quantity 



Q = 2/4/2 



n 
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is a constant and defining X = I^ 2 then implies that the system ([7f)-([T0]) is solved 
by solving the non-linear differential equation 



d 2 X 



dt 



, +oo 2 (t)X 



Q_ 

X 3 ' 



2 



Eq. (jTTj) corresponds to the one-dimensional Ermakov system and is known as 
the Pinney (or Ermakov-Pinney) equation fUJ ■ The general solution of Eq. (jllj) 
is given by [TT] 

X 



AX\ + BX\ + 2CXiX 2 J 1/2 , (12) 
where Xi^it) are linearly independent solutions to the equation 

d 2 Xi 



+ u\t)Xi = 0, (13) 



dt 2 

the constants {A,B,C} satisfy the constraint 

AB ~ C 2 = §2 (14) 

and 

W = Xl —-X 2 — (15) 

is the Wronskian. It follows, therefore, that the evolution of the wavepacket for a 
given trapping frequency, uj 2 (t), is fully determined by solving the linear oscillator 
equation (|T5j) . 

We now proceed to consider cosmological dynamics. Our aim is to identify a class 
of models where the expansion (or contraction) of the universe can be represented 
as an Ermakov system. Recently, it was shown that the Einstein field equations 
for a spatially flat, isotropic and homogeneous Friedmann-Robertson- Walker (FRW) 
universe can be reduced to an Ermakov-Pinney equation when the matter is a mixture 
of (non-interacting) perfect fluids [13j . However, in this case, the sign of the non- 
linear term in the Ermakov-Pinney equation is opposite to that of Eq. (jllj) . This 
difference in sign arises because the energy densities of the matter fields are assumed 
to be positive-definite. 

In view of this, we consider the effects of introducing spatial curvature into the 
analysis. Einstein's field equations for the spatially curved FRW cosmology with 
perfect fluid matter are given by 

1 / A„ \ 2 O U 

(16) 

111 I o u~ 

^ + 3H(p + p)=0, (17) 

where H(r) is the Hubble parameter, a(r) is the scale factor of the universe, r rep- 
resents cosmic (proper) time, k is the curvature constant and units are chosen such 
that Newton's constant is given by AitG = 1. The pressure and energy density of the 
matter are given by p(r) and p(r), respectively. In principle, Eqs. (j!6j) and (j!7j) are 
solved once the equation of state, p = p(p), has been specified. Eq. (j!7j) follows as a 
direct consequence of energy-momentum conservation. 
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Eqs. (fTT)|) and ()17j) may be expressed as an Ermakov system by interpreting the 
matter source as a self-interacting scalar field, <fi = <f)(r), that is minimally coupled 
to Einstein gravity. Formally, this is equivalent to expressing the energy density and 
pressure as 

'=H£) ,+t ™ * =1 2@f-™ (i8) 

where U(4>) represents the potential energy associated with the self-interactions of 
the scalar field. Differentiating Eq. |T6|) with respect to cosmic time and substituting 
Eqs. (fT7|) and (fT%|) then implies that 

1 d 2 a 1 / da\ 2 I ' d<p\ 2 k 



and defining a new independent variable 



d d , . 

— = a— 20 

dr dt y ' 

transforms Eq. (|19j) into the Pinney equation: 

d 2 a I d(j)\ 2 k 

He + \Tt) a= ar (21) 

A direct comparison between Eqs. (|TT|) and (|21|) indicates that the dynamics of a 
positively curved (k > 0) FRW cosmology can be modeled in terms of a harmonically 
trapped Bose-Einstein condensate when cosmic time, r, is related to 'laboratory' 
time, t, through Eq. ([20)1 . The width of the wavepacket plays the same role as that 
of the cosmological scale factor: 

h < — ► a 2 (22) 

and it follows from Eqs. (J7|) and (J2~U|) that the radial momentum of the wavepacket 
is associated with the expansion velocity of the universe: 

da 

/, * * 2-. (23) 



This implies that the Hubble parameter may be identified with the combination 

J/> — if (24) 

and consequently the comoving Hubble radius is determined by the inverse of the 
wavepacket momentum: 

5? ~ !• (25) 
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Cosmological Parameter 


Wavepacket Parameter 


a 


jl/2 


H 2 




P 


3J 4 


V 


I 2 uj 2 - 3/ 4 



Table 1: The dynamical correspondence between a positively curved FRW cosmology with a perfect 
fluid matter source and a harmonically trapped, quasi-two-dimensional Bose-Einstein condensate. 
The width, radial momentum and energy of the wavepacket are analogous to the scale factor of 
the universe, the Hubble parameter and energy density of cosmic matter, respectively. The cosmic 
dynamics is determined once the relationship between the pressure and energy density has been 
specified. This is equivalent in the condensed matter context to choosing the time dependence of 
the trapping potential. 




Furthermore, substituting Eqs. and (jUj) into Eq. (fTUj) and comparing the 

result with the Friedmann equation (jl6|) implies that the energy of the condensate 
wavepacket corresponds to the energy density of cosmic matter 

h <— ► |, (26) 

where we identify Q = k. Finally, the kinetic energy of the scalar field is related to 
the trapping frequency of the condensate such that 

(27) 

and this implies that the pressure of the field is determined by 

p < — ► I 2 uj 2 - 3/ 4 - (28) 

Thus, the width, radial momentum and energy of the condensate wavepacket can 
be identified with a corresponding cosmological parameter. This correspondence is 
summarized in Table 1. Eq. may now be viewed as a novel way of expressing 
the Friedmann equation in terms of properties of the wavefunction. Moreover, by 
employing the correspondences in Table 1, we deduce that Eq. © is equivalent to 
the conservation equation (fT7jl . 

It is of interest to consider specific classes of cosmological models and their con- 
densate analogues. For example, inflationary solutions are of particular importance 
to early universe cosmology. The condition for inflationary expansion is that the scale 
factor of the universe should accelerate (with respect to cosmic time) and it is straight- 
forward to verify that this condition is satisfied whenever the radial momentum of 
the wavepacket grows: 

d 2 a dh , , 

— > — - > 0. 29 
dr 2 dt v ; 
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Another important class of models is characterized by the barotropic equation of 
state, p = (7 — l)p, where < 7 < 2 is a constant. In this case, the correspondence 
in Table 1 implies that the trapping potential satisfies u 2 = 37/4/12 and substituting 
this condition into Eq. (^J) and integrating implies that 



.4 



7-37/2 ' 

1 2 
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37^ 

7(2+37)72' 
1 2 



(30) 



where A is an arbitrary integration constant. A further substitution of Eqs. (JJJ) and 
(J3UJ) into Eq. (fTUj) then yields a first-order differential equation for the width of the 
wavepacket: 

dh _ r 

(61) 



^ = 2^2Alt 3 " )/2 - 1, 
at 



where we have specified Q — 1 without loss of generality. Eq. (j3*Tj) may be integrated 
to yield the solution 




r (3 7 +2)/4 



(3 7 + 2) At, 



(32) 



where 2-^1 is the hypergeometric function [21], At = t — t Q and to represents the second 
integration constant. 

Eqs. (J3U]) and (|3*2^1 together determine the time dependence of the condensate 
trapping potential that mimics a cosmology with a given equation of state, 7. The 
limiting case of 7 = represents a universe dominated by a cosmological constant 
and corresponds to a vanishing trapping potential. On the other hand, a universe 
dominated by relativistic matter has an equation of state given by 7 = 4/3 and in 
this case Eq. (jS^j) simplifies to 



sin X \fx — V x — x 2 



1 



A 

{AAfl^ 1 



At 



2A to 2 / 3 ' 

Eq. ()32j) also simplifies for a universe dominated by pressureless matter (7 = 1): 

1 



(33) 



Ssin- 1 ^ - (3 + 2y)^y - y 2 = —At 

(3A) 1 / 5 1 
2A a; 2 / 5 ' 



y 



(34) 



Although Eqs. (J3l?|) and can not be inverted analytically, the explicit time depen- 
dence of the trapping potential can be expressed in a closed form for a 'stiff' perfect 
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fluid (7 = 2). We find in this case that 



2A 



4(At) 2 
3A 



1/2 



to 



2 [A- 2(At)'< 



l2 ' 



(35) 



To summarize thus far, it has been shown that positively curved, perfect fluid 
FRW cosmologies can be modeled dynamically in terms of quasi-two-dimensional 
Bose-Einstein condensates, where there exists a one-to-one correspondence between 
the type of matter in the universe and the functional form of the time-dependent 
trapping potential of the condensate. The physical properties of the wavefunction 
can be identified with the fundamental cosmological parameters. 

It is of interest to consider extensions of the correspondence to three-dimensional 
condensates. The moment method as outlined above can not be applied in three 
dimensions unless certain approximations are made jTH]- However, an alternative 
approach is to employ the well known scaling properties that are exhibited by the 
wavefunction under the evolution of a time-dependent confining potential [2(Il l2*T| 12*5] . 
In this approach, the ansatz 



1 



Xexp 



i b, 



2 b 



t 2 
r 



(36) 



is made, where b = b(t), x — x(^ r ) an d a subscript denotes differentiation with re- 
spect to laboratory time, t [20] • Substituting Eq. (f3*6*|) into the non-linear Schrodinger 
equation ((TJ and defining new variables 



a = - 



0= Z 



dt 
b 1 



(37) 



then implies that Eq. (JTJ) reduces to 

i dx 
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1 u 2 a 2 x 



b 2 de 



2b 2 

1 



b 3 

ld 2 b 
bdt 2 



b 2 
b 4 



(3? 



where io>o is a constant and the last terms on both sides of Eq. (|38|) have been 
introduced by hand. 

If the interaction strength, g, is time-independent, it is not possible in general to 
separate the left and right hand sides of Eq. (f3*H|) . On the other hand, the scattering 
length may be varied by employing a magnetic field induced Feshbach resonance |26j . 
If the interaction is tuned such that g(t) oc b(t), all terms on the left-hand side of 
Eq. (f3*H)) scale as b~ 2 and, in this case, Eq. (f3*Hjl separates when b(t) satisfies the 
Ermakov-Pinney equation 



d^ + ^ b = ¥- 



(39) 
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It then follows that the rescaled wavefunction, x, formally satisfies the same non- 
linear Schrodinger equation as the original wavefunction but with a constant trapping 
frequency and interaction strength. Comparing Eqs. (|2Tjl and (|39|) implies that the 
scaling parameter, b, may be identified with the scale factor of the universe, where 
laboratory time is again related to cosmic time through Eq. (|2U|). 

We conclude by highlighting a number of questions that would need to be ad- 
dressed in order to realise the above mathematical correspondences in a laboratory 
environment. In general, a sonic horizon may form in an expanding Bose-Einstein 
condensate, where the speed of the fluid exceeds the sound speed. Such an horizon 
forms due to an effective velocity field that increases with distance from the centre of 
the condensate This could be problematic from a dynamical point of view since 
it may introduce physical effects into the system that may not be intrinsic to the 
cosmological model under consideration. In particular, ambiguities may arise when 
interpreting the analogue of cosmological particle production in a time-dependent 
gravitational field. 

A key assumption that was made in establishing the correspondence between 
the condensate and cosmological systems was that the dynamics of the condensate 
wavefunction can be described in terms of the Gross-Pitaevskii equation ([TJ at each 
moment of time, i.e., that the configuration reacts instantaneously to changes in the 
trapping potential and scattering length of the atomic interactions. If this assump- 
tion is to remain valid, the majority of the atoms must remain in the condensate 
state (mean-field approximation) and the particle density and scattering length must 
be sufficiently small (dilute gas approximation). Moreover, to avoid complications 
with delays in the system, the timescales over which external parameters such as the 
trapping potential change should exceed the timescale for a pair of atoms to fully 
interact jHllHZI- Overall, these requirements imply that the trapping potential should 
change at a sufficiently slow rate. In this case, the correspondence (|27|) implies that 
only those cosmologies where the kinetic energy of the scalar field is sufficiently small 
could be simulated 1 . 

A further important scale in the condensate is the 'healing' time, theai- The healing 
time is the characteristic timescale for variations in the wavefunction and plays an 
analogous role to that of the Planck time in cosmology. Comparison of cosmic and 
laboratory time through the correspondence (J2U|) implies that adr = dt > theai and, 
for a cosmic timescale of the order of the Hubble time dr = H^ 1 (representing the 
size of the cosmic horizon), this implies that a > th ea iH. Such considerations suggest 
that the correspondence we have developed could only be realised in practice for 
cosmological models satisfying a lower bound on the scale factor. 

Nevertheless, one advantage of establishing correspondences between cosmology 
and condensed matter physics through Ermakov systems is that insight into the hid- 
den symmetries of the two systems may be uncovered. For example, in the absence 

1 It should be noted, however, that in an inflationary context, the scalar field is necessarily rolling 
very slowly down its interaction potential. 
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of a trapping potential, the action for a two-dimensional Bose-Einstein condensate 
described by the non-linear Schrddinger equation (~~~|) is symmetric under a global 
SL(2, R) reparametrization of the time parameter: t = (at + f3)/(^t + 5), where 
{«,/?, 7, 5} are arbitrary constants satisfying the constraint a5 — (3^ = 1 [23 129*1 - 
In general, this symmetry is broken by the inclusion of a time-dependent trapping 
potential. However, for the specific case where the potential varies as u(t) oc t~ 2 , 
there exists a discrete duality symmetry with a = 5 = and (3 = — 1/7 "2H]. This 
implies that t = —f3 2 /t and such a transformation maps an expanding solution onto 
a contracting one and vice-versa. 

The associated cosmological model for this trapping potential may be deduced in 
the limit where the spatial curvature is dynamically negligible, as is the case during 
inflation. It follows from the correspondence summarized in Table 1 that for this 
trapping potential the scale factor has a power law dependence, a(r) oc r n for some 
constant n, and the scalar field varies as 0(r) oc In r. This solution plays a central role 
in inflationary cosmology because it is one of the few models where the density and 
gravitational wave perturbation spectra can be calculated exactly without invoking 
the slow-roll approximation "2U] (for a review, see, e.g., Ref. jSI]). In effect, the 
solution is a scaling solution since the kinetic energy of the scalar field and the Hubble 
parameter decrease at the same rate. In view of its significance, it is interesting 
that this cosmological model is selected directly from symmetries exhibited by a 
class of Bose-Einstein condensates. It would be interesting to investigate whether 
an associated symmetry exists in the cosmological context as this would provide 
further insight into the density and gravitational wave perturbations generated during 
inflation. 

Finally, an SL(2, R) symmetry also arises in superstring-inspired inflationary cos- 
mologies such as the pre-big bang scenario [22], where a duality symmetry relates 
contracting and expanding solutions. (For reviews, see, e.g., Ref. (22] )■ In these 
models, the dilaton and axion fields in the Neveu-Schwarz/Neveu-Schwarz sector 
of the string effective action parametrize the SL(2, R)/U(l) coset and the existence 
of this symmetry implies that the density perturbations arising from fluctuations in 
these fields can be calculated exactly [M|. Furthermore, inflation in pre-big bang 
cosmology can be re-interpreted as the contracting phase of a universe dominated by 
a stiff perfect fluid [215] and this corresponds in the context of the present work to the 
condensate given in Eq. (""""")) . 
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